Thermal states of the Kitaev honeycomb model: a Bures metric analysis 
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We analyze the Bures metric over the canonical thermal states for the Kitaev honeycomb mode. In this 
way the effects of finite temperature on topological phase transitions can be studied. Different regions in the 
parameter space of the model can be clearly identified in terms of different temperature scaling behavior of the 
Bures metric tensor. Furthermore, we show a simple relation between the metric elements and the crossover 
temperature between the quasi-critical and the quasi-classical regions. These results extend the ones of i29ll30ll 
to finite temperatures. 
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INTRODUCTION 

In the last few years the use of the notion ^t/eZ/fy between 
quantum states to analyze quantum phase transitions (QPTs) 
has proven to be very fruitful [1]-|26]. The motivation comes 
from the fact that a major change in the ground state of a 
many-body system takes place near a quantum phase transi- 
tion [27]. This statement can be made quantitative by bor- 
rowing concepts from quantum information, in particular, dis- 
tinguishability measures between quantum states |28J. The 
simplest of such measures is provided by the amplitude of the 
overlap between the two states. The quantum phase transition 
can be detected by a drop in the fidelity between two ground 
states corresponding to slightly different values of the param- 
eters defining the system. Moreover, this procedure gives rise 
to an associated metric tensor g in the parameter space asso- 
ciated with the system Hamiltonian, whose divergences can 
signal quantum critical points. This strategy has been suc- 
cessfully applied to many different systems. In particular Refs 
dig], [29]-[31] deal with quantum phase transitions involving 
topological order ll32[l . where a local order parameter does not 
exist. 

One can further generalize this fidelity approach to study 
how temperature affects the quantum phase transition. This is 
done by considering the Uhlmann fidelity fss'] between Gibbs 
states corresponding to the many-body quantum system. In 
lfl4l [34I1 it has been shown that different regions in the pa- 
rameter space, called quasiclassical and quantum critical, can 
be clearly identified by the different scaling behaviors of the 
resulting Bures metric tensor g with temperature. 

In this present work we study the effects of temperature 
on the topological phase transition in the Kitaev honeycomb 
model [35]. We extend the results of lEgilsotl to finite temper- 
ature by considering the metric tensor in the parameter space 
that comes from the Uhlmann fidelity between nearby Gibbs 
states. We find indeed that different regions in the parame- 
ter space can be identified according to different scalings in 
temperature of the metric elements, and moreover, we point 
out that the crossover behavior near the critical point can be 
identified by considering a ratio of metric elements. 



Let's begin by reviewing some of the background material 
supporting the idea that near a quantum phase transition the 
degree of distinguishability between mixed states of a many- 
body Hamiltonian is enhanced. The mixed states p we will be 
interested in depend on a set of parameters, e.g. tunable cou- 
plings, {A^} that defines the hamiltonian of the system. The 
problem of distinguishing nearby quantum states can then be 
recast into a problem of estimating the value of the parameters 
{Ap} [36]. The estimation is made from the probability distri- 
bution of the possible results of measurements performed on 
the many-body states p({Ap}). To be concrete, let's assume 
that generalized measurements {E{ei)} are performed on the 
states /3({A^}), with possible results {e^}. This gives rise to 
the probability distribution p(ei I {A^}) = TT{E{ei)p{{X^)}) 
for Ci, given the parameter {A^}. 

A natural distance, called Fisher-Rao distance, can be in- 
troduced in this probability space and is given by 
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which induces a metric 17^1, onto the parameter space {A^} 
given by 



9tLl^ = 
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The Bures distance between two density matrices p and 
p + dp is the natural distinguishability measure that arises as 
the maximization of ([T]l with respect to all the possible gener- 
alized measurements 



dsl 



AP,P + dp) = max{£; 
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Following Refs 1IT4I I34I1 this last expression can be ex- 
panded in the following way: ds^ = J2fiu 9%d,Xf_idXi, + 
g'^^dX^dXv, with the classical and nonclassical metric ele- 
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ments g^^ and g'^l given by 
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where we have used the diagonal basis for p, i.e. p = 
with Pi its eigenvalues, together with the depen- 
dence of these states on the set {A^}. We see that the line 
element separates into a classical and a nonclassical part, the 
first one depending only on the probability distribution {pi}, 
while the second one depends on the set of eigenstates {\i)} 
themselves. The classical part is still the Fisher-Rao distance 
between nearby probability distributions pi and pi + dpi, after 
the maximization over the generalized measurements is per- 
formed JstIi . The expression for these measurement operators 
is known, and turns out to be a projective measurement asso- 
ciated to the Hermitian observable [38] 
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Let's consider the expectation value of M for the thermal 
state p: 



Tr(pAf) - Tr^UpiP + dp)^P = 1 



1 



-dsi 
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after expanding this expression in dp |36]. Since the metric 
elements diverge at the critical point for T = 0, or have a 
finite peak for T 7^ 0, we see that in principle, by measuring 
this observable it would be possible to detect the quantum 
phase transition and ge tting an optimal estimation of the 
Hamiltonian parameters |15j. Notice that M depends on the 
(possibly unknown) state p; therefore the experimental real- 
ization of this measurement may require iterative strategies 
I fSD . We point out that the divergence of the metric elements 
across a QPT is a sufficient, but not a necesary condition ||7(]. 

In the next section we will obtain the explicit form of the 
Bures metric elements for the thermal states of the Kitaev hon- 
eycomb model. 



BURES METRIC AND THE KITAEV HONEYCOMB 
MODEL 



The Kitaev honeycomb model 13511 is one of the most im- 
portant examples of a solvable model in two dimensions that 
exhibits topological phases, with non-abelian anyons. We will 
study the effects of temperature over its critical lines, by con- 
sidering the Bures metric for nearby thermal states. We will 
first briefly review the model and its diagonalization and then 
compute and analyze the associated Bures metric. 



The honeycomb model 

The model consists of spins 1/2 placed at the sites of a hon- 
eycomb lattice, with nearest neighbor interactions. There are 
three types of links at each site, called x, y, or z, depending on 
their direction in the lattice. The Hamiltonian of the model is 
given by 



H = -J,, 
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where the ends of the corresponding links are labeled by j, k, 
and Jj., Jy and are the corresponding model parameters. 
The phase diagram can be separated basically into two regions 
of the parameter space, given by the following inequalities 

\JA < \Jx\ + \Jyl \Jy\ < \JA + \J^l \Jx\ < \Jy\ + \Jzl (8) 

with the equal signs signaling the lines of quantum criticality. 

The region outside these inequalities is gapped and contains 
abelian anyons, while the region inside is gapless and gives 
rise to non-abelian anyons in the presence of a magnetic field. 
To diagonalize this model, we replace the spin 1/2 at every site 
i by four Majorana fermions, by writing each Pauli operator as 
(7° = ibjCj, with b'^ and c Majorana fermions, and a = x,y,z. 
These Majorana operators act on a 4-dimensional Fock space. 
This enlarged space has redundant degrees of freedom which 
can be projected onto a physical 2-dimensional Hilbert space 
C suitable for a spin 1/2 by imposing the condition: \tp) E 
C D\^) — where _D is a projector given by Z? = 

b^b'^b^'c. Within this subspace, the Hamiltonian (|7]i can be 
written as 



(9) 



where the index aji, is x, y, or z depending on the direc- 
tion of the link connecting sites j and fc, and the operators 
Ujk = ib'^j^'' b'^l^'' are such that [ujkTH] = 0, [ujkTUpq] — 
and u^jf. = 1. The entire Hilbert space can be decomposed 
into eigenspaces of the operators ujk, indexed by the eigen- 
values Ujk = ±1. It can be proven that the configuration of 
the Ujk that minimizes the ground state energy is given by the 
vortex-free configuration, that is, Ujk — 1 for all links {j, k). 
Given the translational symmetry of this configuration, one 
can perform a Fourier transform, obtaining the Hamiltonian 
in momentum space within a unit cell 

E *^^7(P)4(P)C7(P): (10) 

where p = [px.Py), ca(p) = ^ Er e"'P''"cA(r) and A(p) 
is a 2x2 matrix given by 



zA(p) = 



i/(p) 
^z/(p)* 



(11) 
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with/(p) = e(p) + iA(p), e(p) = 2{J^cospx + Jy cos py + 
J^) and A(p) = 2{Jj, sinp^; + Jy sinpy). 

In these equations r is the position of the unit cell, while 
A, 7 are the indices for the sites inside each cell. N is the 
number of sites in the lattice, and we choose N = 2L^, with 
L odd. The momenta take the values Px{y) — ^x^i*^ — 



L-l 
2 



L-l 
2 



We can further put the Hamiltonian ( fTOl i in a quasi-free 
fermion form by introducing the fermionic operators [39] 



5(p) = ^(ci(p)+ie"^C2(p)) (12) 
6t(p) = i=(cl(p)-ze-'^4(p)), (13) 



with = arg(/(p)). Then the Hamiltonian takes the follow- 
ing diagonal form: 



i7 = ^A(p)[6t(p)6(p)-- 



(14) 



where A(p) — |/(p)| — y^e(p)M^"A(p)2 represents the 
quasiparticle excitation energies. 

The ground state is given by \gs) — Y\p (ci(p) + 
ie~*^C2(p)) |0), with |0) the vacuum state of the Majorana 
operators ca(p). The excited states are created by applying 
6^(p) onto the vacuum state. 



Bures metric over thermal states 



We are interested in characterizing the effects of finite tem- 
perature on the quantum phase transitions of this model. We 
consider then the Kitaev model in thermal equilibrium with 
a heat bath at temperature T, and obtain the expression for 
the Bures metric (|4]l induced onto the 4-dimensional param- 
eter space of the Kitaev model {/3, J^}, with a = x,y,z 
and (3 = 1/T. In this situation, its states will be described 
byp(T,{J4) = Z-^e-f'" = Z-^j:^ef'^^\i)m, with 
Z = TTe~^^ , Ei and \i) the eigenvalues and eigenvectors of 
the Hamiltonian. Using the equations (|4]l, the expression for 
the fermionic operators (fT3] i. together with the general results 
in lfl4ll . we obtain the following classical metric elements in 
the thermodynamic limit 
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^cosh(A/3) + 1 
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1 



A^ 

cosh(A/3) + 1 
1 

cosh(A/3) + 1 A2 
■ ^^-^^^ cosh(A/?) + l (l)' 

1 ^x^y 

- cosh(A/3) + 1 

1 2e 
cosn(A/?) + 1 A-^ 



(15) 



where a e {x,y}, ilx — 2{cos{px)e{p) + sm{px)A{p)) and 
Uy = 2{cos{py)e{p) + sin(py)A(p)). 

The quantum or nonclassical metric elements are 



nc 

yjaJb 



327r2 



dpxdpy 



cosh(A/3) - 1 e^efc 
cosh(A^) + l A4 



(16) 



wherea,6 e {x,y,z}, Qx = 4(7^ sin(pa;) + Jj^ sin(pj;-pj^)), 
Qy = -4( s\n(px - Pv) - Jz sin(py)) and 9^ = -2A(p). 

In references 1401 kill the authors showed that a related 
quantity, the fidelity susceptibility, could characterize the 
quantum phase transition between these two phases at zero 
temperature, by exhibiting a divergence at the critical lines 
(O. Indeed, similar results can be obtained by taking the 
limit T ^ in ( fTSl ) and ( fTSI l. Furthermore, since in this limit 
the metric elements are inversely proportional to the energy 
gap 10], they exhibit a series of peaks in the gap less region, 
at values of the momenta po for which A = 114 ill . In the 
case of finite temperature the metric elements don't exhibit 
divergences but finite peaks at the critical points, while the 
peaks in the gapless region are smoothed out, since all the 
momenta po are occupied with a probability given by the 
Boltzmann distribution. Figure [T] shows a plot of the metric 
element g^j^, for a finite size of iV = 2L2, L = 101 llsoll . 
for r = and T — 0.01. As can be seen, as soon as the 
temperature is finite, the peaks caused by finite size effects 
are smoothed out. 



Temperature analysis of the metric elements 

In the first part of this section we will characterize two dif- 
ferent regimes in the parameter space Jq}, namely, the 



4 




J 

X 

FIG. 1: (Color online). Plot of g^jJN for a finite size of iV = 
2L^, L — 101, as a function of Jx = 2/3 — Jy. The black and 
red curves correspond to T = and T — 0.01, respectively. As 
soon as temperature is finite, the peaks (due to finite size effects) are 
suppressed. 



quasi-classical and quantum critical, in terms of the different 
scaling behavior of the metric elements with temperature. The 
quasi-classical region is located away from criticality, for the 
gapped phases outside the region (O, and the range of temper- 
atures are such that /3A ^ 1, with A the fermionic gap of the 
system, given by A = 2( — — Jj^) for | | > | | + | |, 
and similarly for the other two gapped regions. The quantum 
critical region is located at the lines of criticaUty that separate 
the gapped and the gapless phases, at finite temperature. For 
this last case, we will analyze the scaling behavior of the met- 
ric elements for A = 0, and approaching the quantum critical 
point by taking T ^ 0. 

Let's proceed first with the scaling of the classical metric 
elements in the quasi-classical region. We will focus mainly 
on the region \Jz\ > | Ja: | + 1 •/y | , but the analysis is very similar 
for the other two gapped phases. Since /3A(p) > /3A ^ 1, we 
can do the following approximation: costj(^A)+i ~ 2e~^-'^^P\ 
The exponential e^'^'^'^P' has a sharp peak centered around 
the minima of A(p), which is given by p = (±7r, ±tt). The 
asymptotic behavior of the metric elements with temperature 
can be obtained by performing a saddle point approximation, 
by expanding A(p) to second order, the rest of the integrand 
up to the first nonzero term around the minima points, and 
extending the limits of integration in p from to cxa. This 
results in the following scaling behavior of the classical metric 
elements 



5^fc(/3A » 1) « r"e-^/^, (17) 

with a ~ 1 for g^^, a = for a — f], b — J^, Jy, Jz and 
a = — 1 for a,b = Jx,Jy,Jz- 

A similar analysis can be performed for the nonclassical 



metric elements, with the result 

<7'-(/3A » 1) « /-(T = 0) + /(J,, Jy, J,)T2e-^/^, 

(18) 

with the same exponent, T^, for all these metric elements, 
with a nonuniversal function f{Jx,Jy,Jz)- 

We focus now on the quantum critical region. We will per- 
form the scaling analysis with temperature inside and along 
the lines of criticality that separate the gapped from the gap- 
less phase, and then take the limit T ^ 0. In this limit, all 
classical metric elements vanish, so we are left to analyze the 
nonclassical ones. 

Inside the gapless region (|8]l, the scaling with tempera- 
ture is dominated by the divergence of A^''^ at its zeroes, 

given by — tt ± arccos ( ""^j^j — ~ ) ^"'l Py — ^ 

arccos f ^2j''j " ) ' around which A sa a (px — ) + /? (pj, — 
Py ), with a and /3 constants. The asymptotic behavior of the 
integrals can be well captured if we shift the limits of integra- 
tion from (— TT, tt) to (0, 2tt), and focus on the region around 
iPx^Py)- The factor ^°^h[A^j+i presents a sharp drop at 
{Px:Py) that gets steeper the smaller the temperature T is, and 
away from this point is equal to 1. For the purposes of ob- 
taining the scaling of the metric elements, it can be very well 
represented by the following piecewise function 

cosh(A/3) - 1 _ r p7t2 for p < T 
cosh(A/3) + 1 ~ \ 1 otherwise, 

where we have taken polar coordinates (p, 9), centered around 
(PxTPy)- Performing a Laurent expansion of ^^f , it can be 
seen that the nonclassical metric elements diverge logarithmi- 
cally inside the gapless region 

5^b'(A = 0,T^0)«ln(T). (19) 

This logarithmic divergence is related to the 2D nature of 
the model, and is different from power law behaviors that have 
been reported for the Ising l(14il and XY chains 1.34,1 at finite 
temperature. 

Along the phase boundaries ([8]l, there is always a direc- 
tion in momentum space along which the quasiparticle energy 
A(p) is no longer Unear but quadratic in momentum. For ex- 
ample, along the path Jx = Jy = ^^2 ' ' ^(P) ~ ^ ''')^ for 
Px — ^Py This quadratic dispersion dominates the scaling 
of the nonclassical metric elements with temperature, so that 
if one restricts the double integrals in (fTSI l to a single integral 
along these paths in momentum space for which the disper- 
sion becomes quadratic, one obtains that 

g2^{A^0,T^0)^T-^/\ (20) 

at the boundaries between the phases, dominating the logarith- 
mic divergence. We have checked this argument numerically 
and confirmed that indeed the nonclassical metric elements 
diverge as T^^/^ for T 0. 
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We conclude that different regimes in the parameter space 
{(3,Ja} can indeed be identified by the distinct scaling in tem- 
perature of the Bures metric elements. We stress again that the 
nonclassical metric elements display a logarithmic divergence 
with temperature. 



QUANTUM VERSUS CLASSICAL CONTRIBUTIONS 

We now remind the reader about the interpretation of the 
Bures metric elements. As we have seen, the Bures distin- 
guishability metric depends on two parts, one classical and 
one nonclassical. The classical metric elements character- 
ize the distinguishability between nearby density states due to 
changes in the mixing probabilities pi of the density matrix 
p, while the nonclassical metric elements characterize how 
the distniguishability is enhanced by the changes in the quan- 
tum states |i) themselves. We can therefore consider the ratio 
between the classical and corresponding nonclassical metric 
elements as a figure of merit to determine the relative im- 
pact on the overall distinguishability due to each classical and 
quantum contribution. Figure |2] shows a contour plot for the 
ratio gj^j^/gj'^j^ in the plane {Jz,T), along the trajectory 
= Jy ^ for 0.48 < Jz < 0.52. In the regions away 

from the quantum critical point — 0.5 and small temper- 
atures and at the quantum critical region above Jz = 0.5 we 
have that the nonclassical metric element has a predominant 
contribution to the Bures line element, while for regions in 
between the classical metric element prevails. 

The contour plot seems to capture the crossover behavior, 
with a crossover at T \ Jz - 0.5]^'', with z ^ ly = 1 IH, 
30l 35]. Indeed, using the expressions for the classical and 
nonclassical metric elements developed in L14<1 . one can prove 
that near a quantum critical point ^ (t) > with A being 
the gap. Since near quantum criticality A \g 



with g the coupling constant driving the transition at gc, we 
have that a contour plot would reveal the crossover at T ~ 

\g-gcr- 




FIG. 2: (Color online). Contour plot of the ratio gj^j^/gj'^j^ as a 
function of Jz and temperature T, along the trajectory Jx = Jy = 
^—J-^. There is an enhancement in discrimination due to the nonclas- 



sical metric element in the quantum critical region as well as away 
from it, for small temperatures. The color scale goes from 0.1 (pur- 
ple), to 0.6 (white). 



can be seen for the Ising and XY models at finite tempera- 
ture. It is tempting to compare this behavior with the usual 
crossover phenomena that separate regions in which the fluc- 
tuations of order parameters are thermal or quantum. Nev- 
ertheless, this connection is blurred for the present model, 
since there is no local order parameter for the Kitaev hon- 
eycomb model. An interesting open question for future in- 
vestigations is to see whether information metric techniques 
similar to those exploited in these paper can help in studying 
quantum phase transitions involving different types of topo- 
logical order see e.g., I,4Z1 . 



CONCLUSIONS 

In this work we have analyzed the canonical thermal states 
of the Kitaev honeycomb model by using the Bures metric ten- 
sor We have shown that metric elements can be used to dis- 
tinguish two regions in the parameter space {T, J^}, namely 
the quasi-classical and quantum critical regions, by different 
temperature scaling of the metric elements. A novel logarith- 
mic divergence with vanishing temperature has been obtained 
for the nonclassical part of the metric elements. 

Furthermore, the ratio between classical and quantum parts 
of metric elements has been used as a figure of merit to com- 
pare the classical and quantum contributions to the overall dis- 
tinguishability of nearby states. It is interesting to point out 
that similar behavior for such ratios between metric elements 
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